Indonesian Journal of Electrical Engineering and Computer Science 
Vol. 25, No. 3, March 2022, pp. 1367~1374 
ISSN: 2502-4752, DOI: 10.1159 1/ijeecs.v25.i3.pp1367-1374 O 1367 


Improving model predictive control's optimization for urban 
traffic 


Ilyas Khelafa, Abdenaceur Baghdad, Mohamed El Hachimi 
Laboratory of Electronics, Energy, Automatics and Data Processing, Department of Electrical Engineering, Faculty of Sciences and 
Techniques of Mohammedia, Hassan II University of Casablanca, Casablanca, Morocco 


Article Info ABSTRACT 


Article history: When it comes to decreasing traffic congestion and enhancing mobility, 

traffic forecasting is critical. However, due to the complicated spatio- 
Received Sep 10, 2021 temporal dynamics of urban transportation networks, which are difficult to 
Revised Dec 25, 2021 describe, this task is tough. Using a model predictive controller (MPC) 
Accepted Jan 17, 2022 provides the control of a traffic network's architecture as well as errors in its 

operations. Based on a real-time simulation, a novel, accurate prediction 

controller for urban traffic was presented in this study to estimate the 
Keywords: number of cars at junctions and their waiting duration. Different 
optimization approaches were employed and evaluated to improve the 
MPC's performance. Simulation results demonstrated that the fmincon was 
very robust and could effectively reduce the number of vehicles in the link, 


Model predictive control 
Optimization algorithms 


S model in comparison with other algorithms This study also includes an in-depth 
Urban traffic control analysis of the characteristics of various prediction horizon sets in an MPC. 
Urban transportation network By increasing the prediction horizon, the amplitude of fluctuation became 
more important, but when Np=4, the fluctuations reduced. 
This is an open access article under the CC BY-SA license. 
© BY SA 
Corresponding Author: 


Ilyas Khelafa 

Laboratory of Electronics, Energy, Automatics and Data Processing 

Department of Electrical Engineering Faculty of Sciences and Techniques of Mohammedia 
Hassan II University of Casablanca 

Casablanca, Morocco 

Email: ilyas.khelafa @etu.fstm.ac.ma 


1. INTRODUCTION 

Communities are quickly expanding and changing as a result of population expansion, economy, 
social, and cultural changes [1], [2]. The number of cars has increased as a result of urbanization, particularly 
in major cities, and issues concerning these vehicles have developed [3]. Despite the fact that building new 
roads every year is highly expensive, traffic is still an issue that creates congestion, especially during rush 
hours [4]. Traffic peaks have become one of the most serious issues in metropolitan areas across the world 
[2], [5], increasing traffic flow and imposing enormous costs on society at all levels. Due to its tight link with 
issues such as air pollution, fuel consumption, and security concerns, the relevance of traffic control has 
grown in recent years. In the form of intelligent transport systems, the smart cities concept offers a potential 
answer to all of these difficulties [3], [5]-[10]. These are usually made up of several connected intersections 
and rely on automated traffic control [1], [11]. 

Until now, traffic management procedures have been based on traditional statistical control 
approaches. However, due to the complicated traffic management system, these approaches are unable to 
provide an acceptable response. As a result, the development of smart approaches can aid in the resolution 
and improvement of urban management issues [8], [12]-[14]. 
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In the planning of transportation infrastructure, traffic modeling and simulations are important [11], 
[15], [16]. Many research projects focus on designing road networks and junctions, evaluating traffic 
circumstances to minimize congestion, decreasing vehicle delays, and increasing road safety. 
Microsimulation traffic models have the potential to greatly change the quality of urban road network 
planning and design [17]-[19]. They may be utilized as part of the decision-making process when comparing 
and analyzing alternative options, when considering changes in traffic regulation, or when optimizing traffic 
conditions for an urban road network in general. The VISSIM microsimulation model was utilized in this 
study [20], [21]. VISSIM is a micro-simulation model for traffic analysis that is stochastic and discrete. 

Different traffic control dynamical models have been discussed. To improve vehicle stopping times, 
traffic capacity, and vehicle delay, Silva et al. developed a multi-objective signal-timing optimization model 
[13]. The model optimized using such approaches increased traffic capacity, according to the testing data. 
The time constant of an intersection control technique was either local or coordinated, and the 
implementation of their method was configured as offline during specific hours of the day for optimization 
[22], [23]. Traffic lights at junctions were regarded one of the most successful ways of regulating urban 
traffic [11], [15], [23]. Travel times and vehicle queues can be reduced along a route with traffic-reduction 
and better-optimized control systems [7], [24], [25]. Fixed-time and smart techniques, can be used to regulate 
traffic [26]. For that, new dynamic multifunctional optimization systems for the predictive control of existing 
network signals have been created and presented as a result of these advancements [7], [27]. 

The purpose of this work is to develop new dynamic multifunctional optimization algorithms for 
predictive traffic control in cities utilizing real-time simulation. This research also aims to show how 
prediction horizon affects model predictive control (MPC). The remainder of the paper is as follows: first, a 
brief explanation of the urban traffic model used during this work is presented. Then, the control algorithm 
and the model predictive control are detailed, with a description of the optimization methods proposed. 
Section 3 presents the simulation results, and the effect of prediction horizon on the MPC. Finally, we 
conclude this paper in the conclusion section. 


2. RESEARCH METHOD 


2.1. Modelization 
The simplified model is presented in this section (called the S model), it is a discrete-time urban 


traffic flow model. The representation’s degree of the real system qualifies the S model, which has a trade-off 
between accuracy and minimal computation. In comparison to previous urban traffic models, the S model 
may incorporate more complicated scenarios such as the time delay required for vehicles running from the 
beginning of the link to the end of the queues, as well as the available storage capacity in downstream 
connections [20]. 

A collection of intersections J and connections L are making up an urban traffic network. A u 
intersection is a place where numerous traffic streams collide. A link (u,d) is therefore a stretch of urban road 
that connects two intersections u and d. The upstream and downstream nodes of a link (u,d) are indicated by 
the letters u and d. For link (u,d), the input and output node sets are u € J and d € O, respectively (e.g., for the 
situation of Figure 1 we have Iua = {iz ,i2, ig } and Oua = {01, 02, 03}). 


a4 01 
| _ | 
= | = | 
=| 21 
-l s 
ER S| 
Se! a 
I =f link (d, u) JI =E link (02, 4) 
. — eee u — _ — _ —_ ee ee eee =] ==? 
io link (i2, u) E = "tt link (u, d) H = itt 02 
E |= 
| P = 
Iž |e 
iz $ 
E 1s 
l 03 


ce) 


Figure 1. Structure of urban traffic network 
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The goal of the S model is to find two states variables at each time step k, nų q and qu a. At each time 
step k, the number of cars in link (u,d) is updated based on the input and output average flow rates over C by: 


Mya (k+1) = Nya (k) - (aka (k) - aua (k)). C (1) 
at every time step k, the number of cars waiting in the line turning to node o is updated over C by: 

quao (k+1) = qu,a,o (K) = (iiao (K) - au,a,o (k)) ~ C (2) 
and at each time step k, the queue length of cars is: 

dua CO = Xoco dudo OO 3) 


We have the following: 

C: Link's cycle time. 

Nua (k): Vehicle’s number during [k.C , (k + 1)C). 

dua (k): Queue lenght during [k.C , (k + 1)C). 

ai a(k): Average link (u,d) entering flow during [k.C , (k + 1)C). 

ala (k): Average link (u,d) leaving flow during [k.C , (k + 1)C). 

a4 a(k): Average incoming flow at the tail of the queue during [k. C , (k + 1)C). 


2.2. Controller design 

Model predictive control is a rolling horizon approach to implement and apply optimum control 
[28]-[31]. An optimum control problem is solved across a prediction horizon in each control step kc, as 
shown in Figue 2, but only the first control sample in the optimal control sequence is executed. After that, the 
horizon is moved one sample, and the optimization is resumed with new measurement data. The optimization 
is re-done using the process's prediction model and an estimate of the disturbance inputs. 

A process model, cost function, and control rule constitute the predictive control structure, as 
Figure 3 demonstrates. By minimizing the cost function and applying the first control signal, the future 
output of a system is predicted over a specified period, using the initial process model of computed control 
signals in the forecast horizon. After that, within the system, calculate it, and repeat the procedure in the 
following phases. The major goal is to anticipate how much time cars spend on the road, which will be 
reduced as future system inputs change at each sample period, therefore improving the system's future 
behavior, which is the green time light. The length of the vehicle queue and the number of vehicles measured 
per link are both inputs to this controller. 
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Figure 2. The concept of MPC Figure 3. The diagram of MPC 


2.3. Optimization 

It is important to choose the efficient optimizer that is used with the nonlinear model in order to be 
able to communicate directly with VISSIM [24], [32], [33]. It's software for simulating multimodal 
transportation operations. In a single model, it depicts all road users and their interactions. All road users are 
realistically modeled using scientifically sound motion models. To examine the reaction of the optimizer for 
this system model, fmincon, nonlinear least-squares solver, particle swarm optimization, and pattern search 
method were used. Matlab is the simulation program that shines out in all methods. The algorithms are 
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explained in the following subsections, along with some suggestions for integrating them into the Matlab 
environment. 


2.3.1. Fmincon 
Finds the restricted nonlinear multivariable function's minimum. The minimum of a problem is 
defined as: 


Ix<x<ux 
y.” SZ 
min f(x) with 4 yeq.x = zeq (4) 
i g(x) <0 
geq(x) =0 


with f(x) is a function that returns a scalar. X, lx, and ux are vectors or matrices that can be passed. y and yeq 
are matrices, z and zeq are vectors, g(x) and geq(x) are functions that return vectors. 


2.3.2. Lsqnonlin 
Solves issues using nonlinear least-squares (nonlinear data-fitting). The nonlinear least-squares 
curve fitting problems is defined as the form: 


min||f (x) If = min (fı (x)? Het fh (5) 


the optimization is obtained with a lower and upper bounds lx and ux on the components of x. x, lx, and ux 
can be vectors or matrices. 


2.3.3. Particle swarm 
The objective of particle swarm optimization is to create a vector x that achieves a local function 
minimum, as the form: 


x = particleswarm (fun, nv, Ix, ux) (6) 


where nv represents the number of design variables, of the function. Particle swarm optimization establishes 
lower and upper bounds on the design variables, x, in order to find a solution within the range: 


Ix <x < ux 


2.3.4. Pattern search 
This method aims to find a minimum of function, using pattern search algorithm. A local minimum 
x, to the function that computes the values of the objective function, is found via pattern search, as: 


x = patternsearch (fun, x0, y, z, xeq, yeq, lx, ux) (7) 


x0 represents the vector specifying an initial point for the pattern search algorithm. 


3. RESULTS AND DISCUSSION 
3.1. Simulation 

The numerical experiments finding are presented in this section. The suggested methodology's 
performance is evaluated in this part using VISSIM, a real-time simulator linked to Matlab through a 
communication interface (COM) [20], with 14 connections and two junctions as shown in Figure 1. Links are 
of the same length: 900 m. All connections are assumed to have three lanes. V free = 50 km/h is the free-flow 
speed, and the average vehicle length is 7 meters. For network roads, the storage capacity C (u,d) is 386 
vehicles. 

In each link, the method has been evaluated at various saturation flow rates: 1600 veb/h for left turn, 
1500 vehb/h for right turn, and 1800 veh/h for through turn. Each road has the same turning rates (0.33 for 
every direction). Furthermore, as Figure 4 demonstrates, all simulations are conducted with variable traffic 
demand and producing congested traffic in the simulator. 
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Figure 4. The inflow demand of vehicles 


The four optimization’s methods are simulated on this platform for the performance comparison: 
Isqnonlin, particleswarm and patternsearch. To explore the four optimization algorithm 


performances, their simulation results are illustrated in Figures 5 and 6. They show the variation of the 
vehicle’s number and their queue length results when the traffic demand shows a variation with 2 peaks, as 
the Figure 4 illustrates. It can be seen that fmincon obtains the best final results. 
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Figure 5. Vehicle’s number 


Figure 6. The queue length 


Observing the results in Figures 5 and 6, we can give the following comments: 
Fmincon achieves the best performance among the four optimization methods. The best and average 
results obtained by fmincon are all smaller than the corresponding results obtained by other algorithms, 
respectively, which demonstrates the strong optimization ability of fmincon in reducing the number of 
vehicles in the link. We can say that fmincon adapts well to different traffic volumes. 
Lsqnonlin performs similar to fmincon, but before the 1st peak of congested demand. Lsqnonlin lost its 
stability with increasing demand. As we can see from Figures 5 and 6 that the average results obtained by 
fmincon and Isqnonlin are close to each other before Ist peak of demand. The worst results of Isqnonlin 
are generally better than that of particlswarm and patternsearch. 
Patternsearch performs similar to particleswarm, but it is more stable. As we can see from Figures 5 and 6 
that the results obtained by the two algorithms patternsearch and particleswarm are close to each other. 
More than that, They present better results face of fmincon and Isqnonlin, but it was limited by the 
increasing of vehicle’s demand (1st peak). Furthermore, they keep the worst results even if the demand 
decreases. Based on the simulation results we can conclude that fmincon is an effective tool for solving 
traffic signal optimization problems and can significantly reduce the number of vehicles, and that makes 
the MPC more effective in front of the other controllers, like fixed time controller [20]. 
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3.2. Prediction horizon effect 

As shown in Figure 2, the prediction horizon, Np, is the number of future check intervals that the 
MPC must anticipate when optimizing its movements to check interval k. Then there's an integer that 
indicates how long the MPC algorithm should take to get the controlled variables to their reference values. 
The prediction horizon Np is essential both in terms of the amount of data fed to the algorithm and in terms 
of the optimization problem's numerical viability [34], [35]. 

Figure 7 presents the results for urban traffic control with different prediction horizons. The most 
noticeable point in Figure 7 is that the number of vehicles is almost dependent of prediction horizon. Note 
that when the prediction horizon increase (from Np=1 to Np=3) that yields to vehicle’s number increase too. 
When Np=4, there is a change about the number of vehicles, with this prediction horizon, we have better 
performance front of Np=3, but the best results still with Np=1. 
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Figure 7. The number of vehicles with different optimization algorithms 


In Figure 8, some fluctuations can be seen. These fluctuations are very important for Np=3. We can 
remark that the controller loses its performances. For Np=4, the MPC takes back the stability, even if it’s 
different from other prediction horizons (Np=1 and Np=2). 

For more illustrations, Figure 9 shows control input, which is the green time length. As we can 
remark, even the output (vehicle’s number) has better results with low prediction horizon, but for control 
input, it’s remarkable that performs poorly in comparison with Np=4 (there is more fluctuations). 
Furthermore, it is necessary to find a good trade-off between the performance of the system and allocated or 
available computational complexity. 
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4. CONCLUSION 

In this paper, four modified optimization algorithms, i.e. fmincon, Isqnonlin, patternsearch and 
particleswarm, were proposed and used to optimize the urban traffic signal. To assess the efficiency of 
different algorithms, a Vissim-Matlab integrated real-time simulation was realized, in which the former 
works as a traffic simulator for evaluating candidate solutions and the latter is responsible for the 
implementation of algorithms. Through this simulation platform, the algorithm optimization performances 
were tested. Simulation results on the traffic network demonstrated that the fmincon was very robust and 
could effectively reduce the number of vehicles in the link. In comparison with other algorithms, fmincon 
achieved the best overall performance in different traffic conditions, like varying the traffic volumes. 

Another goal of this paper is to present the effect of prediction horizon, by studying a MPC 
framework. According to the simulation results, when applying different prediction horizon sets in the 
control of urban traffic, good results are obtained in term of vehicle’s number and the total time spent by 
vehicles. By increasing the prediction horizon, the amplitude of fluctuation became more important, since 
Np=4, the fluctuations reduce, and the controller output is more stable. 
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